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Abstract 

An algorithm is presented in which the Colour-Dipole Cascade Model 
as implemented in the Ariadne program is corrected to match the fixed 
order tree-level matrix elements for e"’“e“ ^ n jets. The result is a full 
parton level generator for e'*'e“ annihilation where the generated states are 
correct on tree-level to fixed order in Og and to all orders with modified 
leading logarithmic (MLLA) accuracy. In this paper, matrix elements are 
used up to second order in Og, but the scheme is applicable also for higher 
orders. A strategy for also including exact virtual corrections to fixed order 
is suggested and the possibility to extend the scheme to hadronic collisions 
is discussed. 


1 Introduction 


Perturbative QCD has been very successful in describing many features of multi¬ 
particle production in high energy collisions. There are, however, several problems 
which have not yet been solved, mostly related to the transition between the per¬ 
turbative and non-perturbative description of the theory. Observables involving a 
few widely separated jets are in principle well described with hxed-order pertur¬ 
bative matrix elements (MEs) for producing a few partons. But to make precision 
comparisons with experiments, it is important to understand the transition of these 
partons to observable hadrons. Our best knowledge of this transition comes from 
hadronization models which describes how multi-parton states are transformed into 
multi-hadron ones. But for these models to work reliably one needs also a descrip¬ 
tion of the collinear and soft partons describing the internal structure of widely 
separated jets and the soft partons between the jets. 


To describe soft and collinear partons it is not feasible to use hxed-order pertur¬ 
bation theory. Not only do the MEs for many-parton states become extremely 
complicated but, since the partons are no longer widely separated, the increase 
in phase space introduces large logarithms which compensates the smallness of Og 
and makes the whole perturbative expansion ill-behaved. To describe the inner 
structure of jets, a more practical approach is to use a parton shower (PS) proce¬ 
dure. Here the large logarithms are resummed to all orders at the expense of only 
keeping the leading logarithmic behaviour of the full matrix elements. 


To get a near complete description of multi-particle production it would be de¬ 
sirable to combine the generation of a few widely separated partons according to 
hxed-order MEs with the evolution of these states according PSs and hnally the 
transition into hadrons using a hadronization model. To do this is, however, highly 
non-trivial and so far there exist no general procedure which is entirely satisfac¬ 
tory. The main problem is that one needs a resolution scale to separate the ME 
generation from the PS one. This scale needs to be small enough to beneht from 
the full ME description, but if it becomes too small the hnal result is spoiled by 
non-physical large logarithms involving the separation scale. 


For some special cases such as e.g. 


e' e 


3 jets 


■p ^ 2 -|- 1 jets 


and 


pp ^ W -|- 1 jet |]9 HT^] , there are working procedures for combining ME and PS. 


Recently there has also been suggested a couple of more general procedures [PT^-jlB 
but none of them can be considered to be the hnal word on the matter. In this 
paper the procedure by Catani and collaborators in |l^ is taken as a starting 


point to create a matching between hxed order ME generators and the colour- 
dipole cascade model (CDM) iHi as implemented in the AriADNE event 
generator. The resulting algorithm is also not complete, but it provides a cleaner 
interface between the ME and PS stages of the generation, and it carefully treats 
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the influence of the available phase space on the Sudakov form factors needed to 
correct the ME generation. 

The outline of this paper is as follows. In section Q the algorithm of Catani et ah is 
presented, reformulated to better fit the CDM. Then, in section the reconstruc¬ 
tion of the resolution scales of a partonic state generated with MEs is discussed, 
followed by the description in section ^ of how the Sudakov form factors are calcu¬ 
lated and a step by step description of the whole algorithm in section In section 
^ a few initial results from the new algorithm for e’''e“ ^ n jets with n < 4 are 
presented. Finally, in section 0, some conclusions are presented together with a 
discussion on how the procedure can be extended to also be used for collisions in¬ 
volving hadrons in the initial state. Also the possibility to improve the description 
to get virtual corrections exact to fixed-order perturbation theory is discussed. 


2 Matrix elements and parton cascades 

To describe e’''e“ annihilation into n jets with fixed order perturbation theory 
one needs to introduce a resolution scale to avoid divergencies in the ME. This 
scale can be defined in many different ways, usually connected to a specific jet 
reconstruction algorithm which can be used both on the ME and on the final- 
state hadrons observed in an experiment. Given such a resolution scale Qq and a 
maximum scale, Q, set by the total center of mass energy, we can write down the 
fraction of n jet events for n < 4, Rn{Q, Qo), given by the second order MEs as 

R2iQ,Qo) = 1 + ®sCo,l(Q) Qo) + ®s^ 0,2(Q; Qo)) 

R'iiQiQo) = ttsC'i,i(Q,Qo)(l-l-ttsC'i,2(Q)Qo))) ( 1 ) 

-R 4 (Q)Qo) = tts^2,2(Q) Qo)) 

where i ?4 receives contributions from qqgg as well as from qqq^q' final states. The 
coefficients Cn,m{Q) Qo) are related to the emission of n partons to O(a^), i.e. Cn,n 
corresponds to resolved tree-level diagrams and Cn,m, with m > n, corresponds to 
the sum of virtual diagrams and unresolved emission diagrams. The problem with 
a small resolution scale comes about because the coefficients contains logarithms 
of Q/Qo which, for small Qo, destroys the ag-expansion. In fact R 3 above becomes 
negative for small enough Qq. 

The parton shower approach is quite different. Here gluons are emitted or split 
into qq pairs in according to simple 1 —> 2 (or, as in the CDM case, 2^3) 
splitting functions in an ordered probabilistic cascade. The emissions are hence 
ordered according to an associated evolution scale so that the first emissions has 
the highest scale and subsequent emissions decreasingly lower scales. The definition 
of the evolution scale, p, is different in different shower models, e.g. invariant mass 
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in PythiA and invariant transverse momentum in AriADNE, but common 
for all models is the occurrence of a Sudakov form factor corresponding to the 
probability of having no emissions between two given scales. If asV{p, z, (p) is the 
splitting probability for emitting a gluon from the initial qq pair, the two-jet rate 
above some scale po becomes 

R 2 {Pm,po) = As^{pni,Po) = exp J dp J dz J d(j)V{p,z,(j))^ , (2) 

where pm is the maximum kinematically allowed scale, and the ^ and (j) integrals 
are over all kinematically allowed values. Performing the integrals and expanding 
the exponential, the two-jet rate can be expressed in the same form as the matrix 
element in eq. (|^) 

R2{pmi Po) = 1 + Q^sC'^l(P,Po) + Po) + • • • (3) 


Similarly the three-jet rate according to a parton shower scenario can be written 


as 


R3{pm,Po) = Os / dp [ dz [ d(t)V{p, Z,(t))As2{pm, p)^S3{p, Po), (4) 

J pQ J J 


i.e. the probability of emitting a gluon at some scale p, times the probability that no 
gluon was emitted above (before) p and the probability that nothing was emitted 
below. The latter is also a Sudakov form factor, symbolically written As^^p^po), 
giving the probability of not having any further emissions from the qgq state above 
the scale po and will in general depends on the full kinematics of the three-parton 
state. Still, it can be written as an exponential of minus an emission probability 
integrated over the allowed phase space. It is therefore easily expanded in powers of 
as and, after integrating an expression similar to the one in eq. (|^ can be obtained: 


R^i^Pmi Po) (P,Po) (It Q^sC'i ^2 ^Pmi Po) T . . .). (5) 


Note that, as in eq. ( 1 . the coefficient Cff can become large and negative, but 
since it is the result of an expansion of an exponential which is always positive and 
below one, the R^ will always be positive and hnite. 


For the four-jet rate there is an additional complication since there are both qqgg- 
and qqq’q’-states, but it can be written on the form 


-R4(Pm,Po) = / dpi [ dZl [ d4>iV{pi, Zi,4>i)As2{Pm, Pi)Ass{pi, P 2 ) X 

J pQ J J 


rPi 


)po 


dp2 / dZ2 / d(t)2 Pgg(p2,^2,02)A5q^^-(p2,Po) (6) 


T T’qq(p2, Z 2 , 02)^5 (p2, Po) 


where is the splitting function for emitting an additional gluon from the qgq 
state and Pqq is the splitting function for splitting the gluon in the qgq state into 
a new quark-antiquark pair. 
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Again, the Sudakovs, and , gives the no-emission probability from 

the qggq and qqq’q’ states respectively and can be written as simple exponentials, 
which can be expanded in as and, after integration, an expression similar to the 
one in eq. (|I|) can be written 

^4(Pm,Po) = asC'0(P,Po)(l + asC'0(pm,Po) + ---)- (V 


The coefficients contains only the leading logarithmic parts[| of the corre¬ 

sponding exact ones, Cn,m- To arrive at a combined parton shower and matrix 
element generator would then mean using the parton shower method but replacing 
the first few coefficients with the exact ones. Here, a more modest approach 
will be taken, where only the coefficients corresponding to the resolved tree- 
level diagrams will be replaced by the exact versions. This amounts to generating 
events according to the n-jet rates 


R2(^Pmy Po) ^S2(^PmiPo)i 


rpm TV/TT? 

Rsipm^po) = Os/ dp Cii (p,n,p)A52(p^,p)As3(p,po), 
-^Po 

pm 


( 8 ) 


-Rqggq(Pm) Po) ~ [ '^Pl [ ^P‘i ^2,2(qggq) (P™-) Pi i P2) ^ 

•1 po on 


' PO 


Asa (Pm, Pi) A53 (pi, P2 ) A5qgg- (P2, Po ), 


.Rqqq'q'(Pm; Po) [ *^Pl [ *^P2 C2^2(qqq'q') (P™.) Pi) Ps) ^ 

J PO J PO 


Asa (Pro) Pi) As3(Pi, P 2 ) A5 (^2, Po), 


where the are the differential full tree-level matrix elements (the 2 ; and (j) 
integrals have been left out for brevity). 

The procedure is then to use a matrix element generator to generate exclusive 
n-parton hnal states, to reweight these states with the Sudakov form factors, and 
then add a parton cascade. 

Note, however, that the p scales are defined in a particular parton shower scenario 
and have no equivalents in the matrix element generation, where the notion of one 
emission being performed before another is not well defined. Instead these scales 
need to be reconstructed from the generated n-parton state. 

In reference this is achieved by using the /c^-algorithm pl|,p5|, the Sudakov 
form factors are then calculated according to an analytical form, and the angular- 
ordered parton shower in HERWIG is added with a special veto procedure. The 
reason for this veto procedure is to avoid double counting, since the reconstructed 


^In fact, most parton shower implementations also contain some non-leading terms corre¬ 
sponding to the modified leading logarithmic approximation (MLLA) |2l| - [^ . In particular, for 
the case of the Ariadne program, = Ci^i and Cq^ = Cop 
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scales do not correspond to the ordering variable in the parton shower. The angnlar 
ordering means that one can have a soft emission at a large angle, followed by a 
harder emission at a smaller angle. 


Here the Colonr Dipole Model implemented in the Ariadne program will be used 
for the parton cascade. The ordering variable is an invariant transverse momentum, 
so the emissions are ordered in scale. In addition, a modihed version of the DICLUS 


27,281 jet clustering algorithm will be used to reconstruct the emission scales given 


as the same invariant transverse momenta. In this way many of the problems in 
ref. [ro are completely avoided. The DICLUS algorithm will not only be used 


to reconstruct the emission scales, but to reconstruct a complete dipole cascade 
history for the n-parton states generated according to the matrix element, i.e. also 
the intermediate states which would occur in a dipole shower producing the same 
state. 


These states will be used to evaluate the Sudakov form factors in exactly the same 
way as they are evaluated in the dipole cascade as explained in section In the 
procedure of Catani et al. the Sudakovs are calculated analytically, which means 
that they do not exactly correspond to a no-emission probability since they do not 
take into account that the available phase space for emission off a quark after it has 
emitted a gluon is changed due to recoils. In the dipole cascade the Sudakovs are 
generated dynamically and takes into account the actual phase space available for 
further emissions. To MLLA accuracy the two approaches are, however, equivalent. 


3 Reconstruction of emissions 


The colour-dipole cascade model |^,|^ describes the emission of a gluon in terms 
of dipole radiation from a colour dipole between two partons. The emissions are 
hence described as two partons going to three, rather than one going to two as 
in conventional parton shower models. This means that gluon coherence is auto¬ 
matically taken into account and that the hrst gluon emission in e’''e“ annihilation 
trivially reproduces the full hrst order matrix element. But there are a couple of 
technical details which are particular to the dipole cascade. All partons are always 
on-shell at each step of the cascade. The conservation of energy and momentum 
is achieved since both emitting partons receives a recoil from the emitted gluon. 
The splitting of a gluon into a qq-pair is also treated as if emitted from one of the 
dipoles connected to the gluon, and the parton in the other end will receive some 

m 


recoil in order to conserve anergy and momentum 


Further more, the scale 


of an emission is dehned in terms of a Lorentz-invariant pj_ of the emitted parton 
with respect to the emitting ones. This is dehned as 

(si2 - (mi -L m2)^)(s23 - (m 2 mg)^) 


p1 = 


5i23 


( 9 ) 
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where parton 2 is the emitted one and Sij and Sijk are the squared invariant masses 
of the two- and three-parton combinations. 


The dipole clustering algorithm (diclus) p7|,^ can be thought of as the inverse of 
the dipole cascade. In each step the combination of three jets which have the small¬ 
est invariant p± are clustered together into two (massless) jets. In the procedure 
presented here, the dipole clustering will be used with a couple of modihcations. 
First of all the information available from the generation of the few-parton state 
will be used. Hence a gluon is only considered to have been emitted from the two 
partons to which it is colour connected]^. So for a colour connected qgig 2 q state 
the only possible clusterings are qgig2 ^ qg2 and gig2q ^ giq. Also, a qq-pair is 
reconstructed into a gluon which is made massless by also considering one of the 
partons connected to the q or q. In this way all partons in the reconstructed states 
are always kept on-shell. A qq pair is only reconstructed to a gluon if they are not 
colour-connected (directly or indirectly). 


When reconstructing a 2^3 emission there is an ambiguity in the directions of 
the emitting partons. If the three partons are transformed into their center of 
mass system, the momenta of the emitting partons are easily obtained from energy 
momentum conservation. It is also clear that the two partons should lie in the 
same plane as the original three. But the angular orientation in this plane is not 
determined. The inverse problem is encountered in the dipole cascade where the 
amount of transverse recoil taken by each of the emitting partons is not given by 
the theory. There are different choices made for each kind of emission. In AriADNE 
the choices are as follows: 


For a gluon emission from a qq dipole, one of the quarks retains its direction 
with a probability proportional to the square of its energy (according to the 
prescription in [pl|] ). 

For gluon emission from a quark-gluon dipole, the gluon retains its direction. 
For gluon emission from a gluon-gluon dipole, the transverse recoil is shared 
among the emitters so that the sum of their squared transverse momenta is 
minimized. 

For a gluon splitting into a qq pair, the spectator parton retains its direction. 


Of these, all but the hrst are completely deterministic and can easily be inverted 
and be used in the reconstruction algorithm. For the hrst, the prescription for 
gluon-gluon dipoles is followed instead^. 

^This requires that the ME generator produces states with a dehnite colour topology. Although 
this information is not physical, it is usually provided by ME generator programs. See discussion 
in ref. 

^This will not influence the results in this paper where only ME will be used and all 
qgq —> qq reconstructions will give the initial qq state where the angular orientation is irrelevant 
for the subsequent cascade. 
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Rather than always selecting the three parton conhguration which has the smallest 
invariant p± to be reconstructed, as is customary in jet algorithms, it is possible to 
reconstruct all possible dipole cascade histories. This is feasible since we are dealing 
with only a handful partons. The procedure will then be to choose randomly 
between these different histories weighted with the corresponding dipole splitting 
probabilities in analogy to the strategy in |^. The splitting probabilities will not 
include a running as as in the normal dipole cascade, since a constant as was used 
in the generation of the parton state. The running of Og will be corrected for at a 
later stage. 


It should be noted that some of the histories may consist of sequences of un-ordered 
emissions which cannot have been produced by the dipole cascade, and these are 
excluded from the histories to choose between. In rare cases it is possible that no 
history can be found which correspond to an ordered sequence of dipole emissions. 
In this case one of the “impossible” histories are chosen, but the reconstructed 
scales are modified so that if the scale of one emission is smaller than the scale of 
the subsequent one, the larger scale is chosen for both reconstructed emissions. 


With this procedure it is now possible to reconstruct a dipole cascade history for 
any n-parton state. All intermediate 2, 3, ..., n —1 parton states are reconstructed 
together with the corresponding emission scales p\i, ..., p\n- 2 - The reconstructed 
scales can then be used to correct the MEs for the running of as by rejecting the 
state with a probability 

1 - ^n”.ya.(pL) (10) 

where the ckso used in the ME generation is taken at the cutoff scale p5_c of 
parton cascade^ to ensure that the probability is positive. 


4 The Sudakov Veto algorithm 

The reconstructed scales and states are also used to calculate the correction for 
the Sudakov form factors. Rather than using the approximate analytic expression 
as a weight, we can use the fact that it corresponds to the no-emission probability 
in a specihc region of phase space. 

Consider a three-parton state generated with the 0{as) ME, where the scale of 
the gluon emission has been reconstructed to p\^. The Sudakov form factor is 
then the probability of there being no emission from the initial qq state before the 
gluon was emitted, i.e. at a scale above p\i, and that there is no emission from 
the qgq state between the scale p\-^ and the cutoff in the ME. By making two trial 

"‘Note that the Qo cut used in the ME generation is not necessarily the smallest possible 
invariant p± scale. 
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emissions with the dipole cascade, one from the reconstructed qq state, starting 
from the maximum scale, and one from the ME-generated qgq state starting from 
and rejecting the whole event if the first was at a scale above p\i or the second 
was inside the ME cutoff, the probability of accepting the event is exactly equal to 
the Sudakov form factor. With this veto procedure the proper phase space region 
is taken into account rather than the approximate limits in the analytic form. 

Special care must also be taken when a partonic state has been generated to the 
highest order in ag used for the ME generation. In that case we want to continue 
generating with the dipole cascade from the last reconstructed scale irrespectively 
of the cutoff in the MEs and we should keep the trial emission from the ME 
generated state and continue the cascade rather than vetoing the whole event. 
This corresponds to changing the four-jet rates in (§) to four-or-more-jets rates, 
e.g. 

fl'Agg+iPmi Po) Q!s [ dpi f dp2 C2^2{gg){Pmi Pli P 2 ) ^ (H) 

JPO JPo 

^82 {.Pmi Pi)^53 (Pl) P2) • 

This cannot easily be achieved in the procedure of Catani et ah, since the angular 
ordering of the parton shower in HERWIG, allows for subsequent harder emissions 
which could result in double counting. 


5 The algorithm step by step 

We now have all the ingredients to present the whole algorithm step by step. We 
assume there is a matrix element generator which can generate complete partonic 
states according to the exact tree-level MEs up to 0{a^). This generator is reg¬ 
ulated by a cutoff Qq to avoid divergencies. In principle this cutoff could be in 
the same invariant p±_ used as evolution variable in the dipole cascade, but to be 
completely general we assume that it is instead e.g. a simple cutoff in invariant 
mass of any two outgoing partons. This is allowed as long as there is no double 
counting (or under counting due to empty phase space regions) as is demonstrated 
below. This matrix element generator should then be combined with the standard 
dipole cascade of AriADNE which has a lower cutoff in the invariant p± given by 
p±c. The constant ttso used in the ME generator is taken at the scale p5_c using the 
same Aqcd as in the dipole cascade. 

The whole procedure will now be as follows: 

1. First the number of partons, n < A^, to be generated is chosen according to 
the integrated jet rates Rn from the tree-level MEs. Note that Rn is larger 
than 1 since we do not include virtual corrections. 


2. Then the momenta of the n partons are generated according to the (!?(q!”“^) 
tree-level ME. Afterwards the invariant p\ of the n partons is checked, and 
if anyone is below the state is rejected and the procedure is restarted at 
step 1. 

3. Now, all the intermediate states 52, ..., Sn-i and scales ..., p]_n -2 
corresponding to a sequence of dipole emissions are reconstructed according 
to the algorithm in section 

4. The generated event is rejected and we restart at step 1 with a probability 
given by equation ([T0|). 

5. We now make a trial emission with the dipole cascade from the state S 2 
starting from the maximum scale limited by the center of mass energy. If 
this emission is at a scale above fhe event is rejected and we restart 
from step 1. If not, a trial emission is performed from the state S^, with a 
maximum scale of If this emission is at a scale above p5_2 event is 
rejected and we restart from step 1. This procedure is repeated for all states 
down to Sn-i- If no rejection has been made, a trial emission is made from 
the ME-generated n-parton state starting from the scale p]_n- 2 - There are 
now two cases 

• If n = the trial emission is always kept and the dipole cascade is 
allowed to continue down to the cutoff and the event is accepted. 

• If n < A^, and all parton pairs pass the ME cut, Qo, the event is rejected 
and we restart from step 1. If any of the partons fail the cut, the trial 
emission is accepted and the dipole cascade is allowed to continue down 
to the cutoff p\^ and the event is accepted. 

For a pictorial description of the procedure, hgure |I] shows the regions in a symbolic 
two-dimensional phase space, in which the Sudakov veto algorithm is used in the 
case of A^ = 4 and n = 2, 3,4. It is clear that if there is only one emission inside the 
ME cut, it is handled by the ME generator and this region is never populated by 
the dipole cascade. If there are two emissions inside the ME cut, the two hardest 
ones are always handled by the ME generator while additional emissions are given 
by the dipole cascade. In this way there is no double counting between the ME 
generator and the dipole cascade. Note that if the trial emission in e.g the n = 2 
case is outside the ME cut and therefore is accepted as shown in hgure ^ further 
emissions inside the ME cut are allowed without any risk of double counting. 

Although the procedure is to add a dipole cascade to an A^-parton matrix ele¬ 
ment generator, the result is that all hnal multi-parton states are distrubuted as 
if generated by the dipole cascade, except that if the n — 2 (with n < N) hardest 
emissions are inside the matrix element cut, their distribution is described by the 
exact tree-level matrix element. 

More precisely, we have the following possible situations: 
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ME cut Qq 



/ 


Dipole cascade cutoff p 


Kinematical 


^ 

phase space 
boundary 


P±1 





n = 


2 


71 = 


3 


71 = 4 


Figure 1: The integration regions (shaded) for the Sudakov form factors in a sym¬ 
bolic two-dimensional phase space, giving the probability to accept a ME-generated 
2, 3 and 4 parton states in the case N = 4. The evolution variable in the cascade 
is assumed to be along the vertical axis, while the horizontal axis is some rapidity 
or energy splitting variable. The cutoff in the ME is along the diagonal line. 


• If the hardest emission is outside the ME cutoff, all emissions are given by 
the dipole splitting functions, even if softer emissions are inside. 

• If the hardest emission is inside the ME cutoff, but the second hardest is 
outside, the hardest emission will be given by the ME, while all subsequent 
ones are given by the dipole splitting functions. 

• If the two hardest emissions are inside the ME cutoff they are given by the 
ME, while all subsequent emissions are given by the dipole splitting functions 
(irrespective of whether they are inside the cutoff or not). 

In all cases the Sudakov form factors are the same irrespective of whether the 
emissions are given by the ME or the dipole splitting functions. 

It is clear that there should only be a small dependence on the Qq since the only 
change when going outside the cut is that the emissions are governed by the leading 
logarithmic expressions rather than the exact ME and these should be very similar 
for a small enough cut. It is important that the extra p± cut used in step 2 is the 
same for the ME-generated state as for the subsequent dipole cascade, otherwise 
the ME may populate phase space outside of the reach of the cascade. This is, 
in fact, not quite trivially achieved in the cascade itself. The reason is that the 
recoils from an emission in one dipole may push the invariant p± of a parton in a 
neighboring dipole below the cutoff pj_c. To be completely consistent, the dipole 
cascade is therefore used with an extra cut vetoing emissions which push other 
partons below the cut. 
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Figure 2: Possible sequences of dipole cascade emissions from n = 2 and n = 3 
states generated with MEs. 

6 Results 


To check that the procedure works, the simplest thing is to test it for N = 3 since 
the Ariadne dipole cascade is already corrected to match the 0{as) ME. And, 
indeed the new procedure completely agrees with standard Ariadne in this case. 
A less trivial test is to look at the A^ = 4 case. For this a modihcation of the 0{ag) 
ME generator implemented in Pythia is used, stripped down so that it only uses 
the tree-level MEs with a simple invariant mass cutoff. This is then used together 
with the new dipole cascade interface algorithm implemented in Ariadne. 


First we look at some standard event shapes. These are known to be very well 
described by the Ariadne program (see eg. |^). We will not directly compare 


to data since this would not be meaningful without a retuning of the hadroniza- 
tion parameters. Instead we will compare the new procedure with the standard 
Ariadne program on parton level. Figure shows the dependence on the ME 
cut Qo of some standard event shapes which should be sensitive to effects. 

The differences w.r.t. Ariadne are very small and rather than showing the event 
shapes them selves, hgure ^ shows the ratio to the AriADNE results. The fact that 
the differences are so small is not a surprise. Already in it was shown that the 
dipole cascade agrees very well with 0(ag) MEs in most regions of phase space. 
For large Qo fhe new procedure reduces to the standard Ariadne program as ex¬ 
pected. The dependence on the Qq is small but not zero. The largest difference is 
found i the D4 distribution which measures the smallest distance, according to the 
Durham jet-algorithm 


, 25|, between two jets when an event has been clustered 


to four jets. It is a bit worrying that the dependence does not seem to disappear 
even for very small yo = Qq/EI^. This may be expected for the D4 distribution 
when the distance between jets are smaller than but it is also present in the 
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(a) Objateness normalized to Ariadne (b) D4 normalized to Ariadne 






Figure 3: Ratios of parton level event shapes at Eqm = 91 GeV for the new ME 
matching algorithm using different values of Qo w.r.t. to the standard Ariadne 
program. The distributions are (a) oblateness, (b) f-jet resolution value for the 
Durham jet algorithm and (c) p^out- (^) same as (c) but with massless 

guarks and no secondary guarks both for the new ME matching and for standard 
Ariadne. In all cases, the full line is with yo = Ql/Q^ = 0.05, long-dashed: 
yo = 0.02, dashed: yo = 0.01 and dotted: yo = 0.005. 


P±out distribution even for fairly large values. The reason can be traced to the 
treatment of quark masses which is not exact in Ariadne (nor is it exact in the 
tree-level matrix elements in Pythia). In addition there are some uncertainties in 
the treatment of secondary quarks in Ariadne[|. As seen in in hgure P, if all quark 
masses are put to zero and secondary quarks are taken away, the dependency on 
I/O becomes much smaller, but it still does not go away for really small values. A 
possible explanation may be that the dipole clustering routine has some problems 
for very small resolution scales as described in |^, but the dependence is in any 
case much smaller than the uncertainties due to hadronization parameters and the 


®An investigation of the effects of quark masses and secondary quarks will be presented in a 
future publication 34 
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Figure 4: The distribution in the Bengtsson-Zerwas angle on parton level using the 
JADE jet algorithm j^/ with i/cnt = 0.03. In (a) the full line is standard Ariadne, 
the dashed line is the new ME matching algorithm with i/o = 0.01 and the dotted 
line is the tree-level 0(ag) ME-only generator in PythiA with i/o = 0.01. In (b) the 
full, long dashed, dashed and dotted lines are all the new ME matching algorithm 
with hq = 0.05, 0.02, 0.01 and 0.005 respectively. 
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basic parameters in the dipole cascade, p±c and Aqcd- 

To really see the influence of the ME matching one must look at details in the 
correlations between jets. Here we will look at Bengtsson-Zerwas angle [0, i.e. 


the angle between the plane spanned by the two largest jets and the plane spanned 
by the two smallest jets in a four-jet event, which is not at all described by the 
standard Ariadne program. In figure ^ we see that the new matching procedure 
is closer than standard Ariadne to the result from the pure 0{ag) ME generator 
in Pythia. It does not, and should not, exactly reproduce the pure ME approach 
since the correlation is smeared by the subsequent soft radiation. 


7 Conclusions and outlook 

The procedure presented here to correct the dipole cascade in Ariadne to match 
the exact 0{a^) matrix elements works. The only additional parameter needed is 
the cutoff Qq in the ME generation, but the results have been shown to be fairly 
insensitive to the value chosen as long as it is reasonably small. Although the results 
presented here has only been for iV = 4, the procedure is completely general and 
can be applied to basically any iV-parton ME generator. One advantage is that 
the procedure is practically non-intrusive with respect to the ME generator used. 
On the other hand the vetoing technique makes the generation very inefficient for 
very small Qq. 
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There are a few uncertainties in this algorithm, mainly to do with the reconstruc¬ 
tions of the emissions. It is not quite clear what to do with reconstructed emissions 
which are not ordered. The procedure presented in section ^ is not unique and, 
in fact, several other options have been tested. The results are, however, rather 
stable w.r.t. such variations. 

In any case the procedure will give multi parton final states where the first N — 2 
emissions are correct to exact 0{a^~‘^) accuracy and all others are correct to MLLA 
accuracy. 

It could be possible to also include the exact virtual corrections to This 

would involve correcting the tree-level matrix elements with the difference between 
the exact and the leading logarithmic parts of the virtual plus unresolved terms, 
for the iV = 4 case effectively turning R 2 and R 3 in equation (||) into 


R2 {.Pm.) Po) 
Rz iPmi Po) 


As 2 (Pm, Po) + as6C^^^^{pm, Po) + {pm, Po), 

rpm A/fr* 

tts / dpc^f(pm,p)x (12) 

Jpo 

(13) 


where the difference between the full coefficients from the summed virtual and 
unresolved diagrams in the matrix element and the corresponding parton shower 
ones, from the expansion of the Sudakov form factors and the running of ag with 
transverse momenta in the parton cascade, are symbolically written = 

These coefficients are free from singularities, but the cutoff in the 
matrix element generation must be chosen so that the differential n-jet rates will 
come out positive in the whole phase space. Further investigations is needed to see 
if this at all can be accomplished. 


It should be noted that a procedure for combining next-to-leading order matrix 
element with parton showers has recently been proposed by Frixione and Webber 
T7|. However, such a generator will produce negatively weighted events and is 


therefore difficult to use as a general purpose event generator. 


Finally it should be noted that it is possible to use the procedure presented in 
this paper also for collisions with incoming hadrons. The additional complication 
is that the reconstruction of the intermediate states and scales must take into ac¬ 
count that there may both be initial- and final-state emissions, where the former 
depends on the parton density distributions of the incoming hadrons. But other¬ 
wise the procedure would be the same: Use a tree-level matrix element generator to 
generate lowest order and subsequent higher order partonic states, reconstruct the 
possible emission histories according to a parton cascade scheme together with the 
corresponding emission scales, calculate the running of Og, calculate the Sudakov 
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form factors with the veto algorithm using trial emissions from the reconstructed 
states, and continue the cascade. 
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